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Abstract 

We evaluate the mass function of virialized halos, by using Press & Schechter (PS) and/or Steth 
& Tormen (ST) expressions, for cosmologies where Dark Energy (DE) is due to a scalar self- 
interacting field, coupled with Dark Matter (DM). We keep to coupled DE (cDE) models known to 
fit linear observables. To implement the PS-ST approach, we start from reviewing and extending 
the results of a previous work on the growth of a spherical top-hat fluctuation in cDE models, 
confirming their most intriguing astrophysical feature, i.e. a significant baryon-DM segregation, 
occurring well before the onset of any hydrodynamical effect. Accordingly, the predicted mass 
function depends on how halo masses are measured. For any option, however, the coupling causes 
a distortion of the mass function, still at z = 0. Furthermore, the ^-dependence of cDE mass 
functions is mostly displaced, in respect to ACDM, in the opposite way of uncoupled dynamical 
DE. This is an aspect of the basic underlying result, that even a little DM-DE coupling induces 
relevant modifications in the non-linear evolution. Therefore, without causing great shifts in linear 
astrophysical observables, the DM-baryon segregation induced by the coupling can have an impact 
on a number of cosmological problems, e.g., galaxy satellite abundance, spiral disk formation, 
apparent baryon shortage, entropy input in clusters, etc.. 

PACS numbers: 98.80.-k, 98.65.-r 
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INTRODUCTION 



A first evidence of Dark Energy (DE) came from the Hubble diagram of SNIa, showing an 
accelerated cosmic expansion, but a flat cosmology with Q m ~ 0.25, h ~ 0.73 and ~ 0.042 
is now required by CMB and LSS observations and this implies that the gap between Q m 
and unity is to be filled by a smooth non-particle component, whose nature is one of the 
main puzzles of cosmology. (fi m> &: matter, baryon density parameters; h: Hubble parameter 
in units of 100 km/s/Mpc; CMB: cosmic microwave background; LSS: large scale structure.) 

If DE is a false vacuum, its pressure/ density ratio w = Pde/pde is strictly -1. This option, 
however, implies a severe fine tuning at the end of the electroweak transition. Otherwise 
DE could be a scalar field <f), self-interacting through a potential V{4>) , so that 

PDE = Pk,DE + Pp,DE = (ft 2 /2d 2 +V((j)), PDE = Pk,DE - Pp,DE , (1) 

provided that dynamical equations yield pk,DE/V <C 1/2, so that — l/3>w>- 1. Here 

ds 2 = a 2 (r)(-dr 2 + dxdx 1 ) , (z = 1, .., 3) (2) 

is the background metric and dots indicate differentiation with respect to r (conformal time). 
This kind of DE is dubbed dynamical DE (dDE) or quintessence; the w ratio then exhibits a 
time dependence set by the shape of V{<p). Much work has been done on dDE (see, e.g., [3] 
and references therein), also aiming at restricting the range of acceptable io(t)'s, so gaining 
an observational insight onto the physics responsible for the potential V(4>). 

As a matter of fact, the dark cosmic components are one of the most compelling ev- 
idences of physics beyond the standard model of elementary interactions and, while lab 
experiments safely exclude non-gravitational baryon-DE interactions, DM-DE interactions 
are constrained just by cosmological observations. In turn, DM-DE interactions could ease 
the cosmic coincidence problem Q], i.e. that DM and DE densities, differing by orders of 
magnitude since ever, approach one another at today's eve. 



In a number o 
were discussed (4, 



3ap_ers, constraints on coupling, coming from CMB and LSS observations, 
This note aims at formulating predictions on the mass function of 
bound systems in coupled DE (cDE) cosmologies, 
known to improve the original Press & Schechter 



>y using Sheth & Tormen [9J expressions, 
approach. Our final scope amounts to 



finding stronger constraints on DM-DE coupling, arising from a comparison of observational 
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data with our predictions, opening a basic window on the nature and origin of these very 
components. Our analysis will be restricted to SUGRA and RP (Ratra-Peebles) 3] 
potentials 

SUGRA V(4>) = (A a+ 70 a ) exp(47T0 2 /mJ) 

RP V{<t>) = (A Q+4 /0 a ) (3) 

(m p = G~ x l 2 : Planck mass), admitting tracker solutions. This will however enable us to 
focus on precise peculiarities, not caused by the shape of V{4>) but by the coupling itself. 
Let us also remind that, once the DE density parameter Qde is assigned, either a or the 
energy scale A, in the potentials can still be freely chosen. In this paper we show results 
for A = 10 2 GeV; minor quantitative shifts occur when varying log 10 (A/GeV) in the 1-4 
range. The RP potential will be mostly considered to test the effects of varying DE nature. 
The effects of coupling can be seen in the background equations for DE and DM, reading 



+ 2(d/a)0 + a 2 V^ = ^16vrG/3/3a 2 p c , p c + 3(d/a)p c = -jl&irG/Sppcjt ; (4) 

here j3 sets the coupling strength and, all through this paper, we take it constant (and, in 
particular, independent of 0; a different case, which can be physically significant and may 
deserve a separate treatment [t|) with values f3 = 0.05 or 0.20 . In previous work this was 
considered a small strength. CMB data set a limit 3 <^ 0.6; a stronger constraint found 
by 0, by studying cluster profiles in n-body simulations, applies to RP potentials only. 

On the contrary, even for the smaller coupling we considered, we find significant devia- 
tions from uncoupled models in non-linear observables. Such deviations exhibit intriguing 
features, suggesting possible ways out from a number of astrophysical problems, while the 
small coupling strength scarcely affects linear observables. 

As far as the other model parameters are concerned, we take density parameters Q m = 
0.25, Vl b = 0.042, Vl c = 0.208 and h = 0.73. This choice approaches the best fit of a ACDM 
model to available data. A best fit to cDE models will certainly yield slightly different 
parameters, but this implies just minor shifts in our quantitative findings. 

BARYON AND DARK MATTER DYNAMICS 

The essential novel feature induced by DM-DE coupling in non-linear structures is 



baryon-DM segregation. This was shown in 



ll| (paper I, hereafter). The reason why segre- 



gation occurs can be easily illustrated by considering, aside of the eq. (J3J), ruling background 
dynamics, the couple of equations telling us how the baryon (DM) density fluctuations 5 b 
(5 C ) and velocity fields 

a - k ' v c,b 

°c,b = t n ■ (5) 
depend on r (Ti, = a/a; k: wavenumber of the fluctuation considered). They read 

5 C " = -5 C '(1 + % - 2/3X) + + t;P 2 W c + jW, 
rt 16 I 

5 b " = -V(l + %) + jWc + n b 5 b ), (6) 

rt I 

e c ' = -e c (i + %- 2[3X) - hi + \(3 2 )n c 5 c - ;W 6 , , 

it 16 I 

e b ' = -e b (i + ~(n c 5c + n b 5 b ) . (7) 

rt I 

Here ' yields differentiation with respect to a = In a and 

X = y/4^/3<j>/(m p H) . (8) 

Both eqs. (UJ) and eqs. ©-(0) show that, as soon as (3 ^ 0, baryons and DM have differ- 
ent dynamics. Neglecting radiation and any hydrodynamical effects, baryons move along 
geodesies. On the contrary, DM particles feel also non-gravitational forces and, as it hap- 
pens in the presence of any non-gravitational interaction, do not follow geodesies. This is 
why baryon-DM segregation occurs. 

Starting from eqs. the evolution of linear fluctuations in cDE cosmologies was 

studied by [4J. Here we shall put ourselves into a different physical context, by lifting 
the restriction 5 Ctb <C 1, but considering just cases when a full relativistic treatment is 
unessential. Then, dealing with scales well below horizon and with non-relativistic particles, 
the long-range force carried by the DE field <fi can be described through corrections to 
newtonian gravity. They amount to assuming: 

(i) DM particle masses to vary, so that 

M c (t) = M c ( n ) exp[-C(0 - <Pi)}. (9) 

(ii) The gravitational constant between DM particles to become G* = jG. 
Here 

C = ^JiQttG/3 f3 , 7 = l + 4/5 2 /3 (10) 
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A proof is given in Appendix A, following 



121 and 



A newtonian treatment is suitable to study the growth of spherical top-hat fluctuations 
in the DM and baryon components and to show how their segregation occurs in the relevant 
physical cases. In fact, non-linearity starts well after matter-radiation decoupling and the 
top scales to be considered, galaxy cluster scales, lay well below the horizon. 

The evolution of fluctuations in this regime was debated in paper I. Here we shall use 
its results to understand the shaping of mass functions. We therefore start from reviewing 
them and complementing them with new quantitative outputs. 



SPHERICAL TOP HAT EVOLUTION 



Although in all cosmologies, apart sCDM, the system of equations ruling the spherical 
growth requires a numerical solution, in cDE models the numerical approach is far more 
essential. In fact, in all cases apart cDE, the only variable describing a top-hat fluctuation is 
its radius R. It is assumed to expand, initially, at the same rate of the scale factor a, obeying 
an equation similar to the scale factor in a closed model. The greater density inside the top 
hat slows then down the increase rates of R in respect to a, so that the inner density p(< R), 
although decreasing, becomes greater and greater than the average density p. Eventually, 
at a time t ta , when the density contrast A = p(< R)/p attains a suitable value Xi R stops 
and starts decreasing. The equation it obeys would then cause R to vanish within a finite 
time t c . In sCDM cosmologies, x — (Sti/A) 2 and t c /t ta = 2. In other non-cDE cosmologies, 
X and the ratio t c /t ta take just slightly different values. 

After tta, however, the very density inside R is increasing, so that, in any realistic fluc- 
tuation, unless radiation disposes the heat produced by the p dV work, virial equilibrium is 
attained and the collapse essentially stops when the fluctuation radius is R v . In a sCDM 
cosmology this occurs for R v /R ta = 1/2. Taking onto account the symultaneous growth of 
a, the virial density contrast is then = 32% ~ 180. In other cosmologies the R v /R ta 
ratio is slightly different from 1/2 and A„ is mostly above 100, but, in most cases, signifi- 
cantly below 180. A usual assumption is that the time when the system is stabilized into 
a virialized configuration coincides with t c . Let us also remind that while, in sCDM, A„ is 
independent from t c , this is no longer true in all models where the density ratio between DE 
and material components depends on time. 
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Quantitative results for the spherically symmetric growth in ACDM models were given by 
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151 ] - In these cosmologies, 



[13^ ; their extension to models with uncoupled dDE is due to 
the simple system of equations yielding spherical evolution requires a numerical solution. 
Let us now consider the same problem in the cDE case. This problem had been considered 



in Q, 



but its approximated treatment did not allow the authors to focus the main physical 



effects. Following 



are already different 



1 11 , le t us outline, first of all, that linear fluctuations in baryons and DM 
12| and have a ratio 

8b _ , 3fi c 



8 C 3 7 fi c + 4/3X// 

Here X is given in eq. (JBJ), 7 is given in eq. (fTUJl : fi = (S c> b/S c> b)/7i. When considering a 
spherical top-hat fluctuation, whose initial radius RrH,i = Rb = Rc, eq. (fTTjl sets the initial 
ratio between DM and baryon fluctuation amplitudes. 

Initial conditions are set so that both Rb and R c initially grow as the scale factor. However, 
because of the different interaction strength, as soon as non-linear effects become significant, 
Rb starts exceeding R c . Hence, a part of the baryons initially within Rth,% leak out from R c , 
so that DM does not feel the gravity of all baryons, while baryons above R c feel the gravity 
also of initially unperturbed DM layers. As a consequence, above R c , the baryon component 
deviates from a top-hat geometry, while a secondary perturbation in DM arises also above 
R c itself. 

After reaching maximum expansion, contraction will also start at different times for 
different components and layers. Then, inner layers will approach virialization before outer 
layers, whose later fall-out shall however perturb their virial equilibrium. This already 
outlines that the onset of virial equilibrium is a complex process. 

Furthermore, when the external baryon layers fall-out onto the virialized core, richer of 
DM, they are accompanied by DM materials originally outside the top-hat, perturbed by 
baryon over-expansion. 

The time when the greatest amount of materials, originally belonging to the fluctuation, 
are in virial equilibrium occurs when the DM top-hat has virialized, together with the 
baryon fraction still below R c . A large deal of baryonic materials are then still falling out. 
But, when they will accrete onto the DM-richer core, they will not be alone, carrying with 
them originally alien materials. There will be no discontinuity in the fall-out process when 
all original baryons are back. The infall of outer materials just attains then a steady rate. 



In order to follow the dynamical evolution of a systems where each layer feels a different 
and substance-dependent force, a set of concentric shells, granting a sufficient radial reso- 
lution, needs to be considered. In the next section we shall provide the equations of motion 
for Lagrangian shell radii and review the whole expansion and recontraction dynamics. 

TIME EVOLUTION OF CONCENTRIC SHELLS 
The spherical growth 

Dynamical equations can be written by using comoving radii b n = R b /a (c n = R™/a) for 
the n-th baryon (DM) shell. Let M; c be the masses of the n-th layer, M£ c be the masses in 
the same layer in the absence of perturbations, and 5M£ C = M£ c — M£ c be the mass excesses 
for the n-th layer. Baryon shells keep a constant mass in the whole process, while 

M c » = M 6 "(r m )exp{-C[0(r) - <j>(r in ))} , (12) 

according to eq. (JOJ). Similarly, let AM& iC (< r) = M^ c {< r) — Mb. c (< r) be the mass excess 
within the physical radius R = ar and let then be 

AM t n = AM b (< bn) + AM C (< b n ) , AM c n = AM b (< c n ) + 7 AM C (< c n ) , (13) 

including a correction aimed to take into account, besides of the mass variation (eq. fT2*|) . 
also that G* 7^ G in DM-DM interactions. 

A systematic discrepancy will arise between b n and c n and will require taking into account 
only the fraction of the last shell below the radius considered. In paper I it is shown that 
the shell dynamics can then be described through the equations 

c n = - - -C(j))c n - G f- , b n = — b n - G—f - , (14) 

\a ) ac n a ab n 

to be integrated together with the first of eqs. (jlj) and the Friedman equation. 

Until b n ~ c n , it is AM fc ™ < AM", for any j3 > 0. Hence, the gravitational push felt by 
DM layers is stronger. The extra term C(j)c n adds to this push. In fact, the comoving variable 
c n has however a negative derivative (see the lower panels in Figure [TJ. As a consequence, 
c n decreases more rapidly than b n , so that the n-th baryon radius may gradually exceed the 
(n + l)-th DM shell, etc. . As a consequence, the sign of AM; — AM; can even invert and 
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FIG. 1: Evolution of a sample of baryon and DM layer radii, extending up to the top-hat radius; 
upper (lower) panels describe physical (comoving) radii. In all plots the leaking out of the upper 
baryon layer is clearly visible. 

one can evaluate for which value of <fi — (pi this occurs. Once 0(r) is known, also the time 
when this occurs can be found. Until then, however, DM fluctuations expand more slowly 
and mostly reach their turn-around point earlier, while baryons gradually leak out from the 
fluctuation bulk. 

The whole behavior is visible in the Figure ^ for samples of R c ^ and c-b (the ordinate 
label x stands for either c or b). The greatest radii shown are the baryon top-hat radii. 
Solid (dotted) lines yield the a dependence for DM (baryon) shell radii. The actual number 
of radii B% c used in the equations depend on the precision wanted; radii do not need to be 
equi-spaced. To our aims (a 0.1 % precision), ~ 1000 radii were sufficient. If more precision 
is needed, it can be easily achieved at the expenses of using a greater computer time. With 
a single processor 1.5 Gflops PC, our optimized program takes ~ 60 minutes to run a single 
case, with ~ 1000 radii. 
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FIG. 2: Density profiles at different a values for a cDE model with (3 = 0.05. Solid (dotted) lines 
refer to DM (baryons). Notice the progressive deformation of the baryon profile (dotted lines) in 
respect of the DM profile (solid line). 
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FIG. 3: Density profiles at different a values for (3 = 0.2 model. The deformation of the baryon 
profile (dotted lines) in respect of the DM profile (solid line) is stronger than in the previous Figure. 
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The full collapse shown in Figure [T] is not expected in any physical case, as the R n 
decrease stops when virialization is attained. Figs. El El therefore assume that the spherical 
growth stops when all DM originally in the top-hat, and the baryons kept inside it, virialize. 
This figure describes the gradual deformation of the top-hat profile. Already at a = 0.2, 
well before the turn-around, the slope of the top-hat boundary, for baryons, is no longer 
vertical. At a = 0.6 (approximately turn-around), not only the baryon boundary is bent, 
but a similar effect is visible also for DM. The effect is even more pronounced at cw ~ 0.92. 
Figures similar to the upper plots in Fig. [T] and to Figs. El El were already shown in paper I, 
although for different model parameters. 

The increased density of DM shells outside the top-hat is relevant, here, because it 
fastens the recollapse of outer baryons. In turn, the enhanced baryon density, outside from 
DM top-hat, modifies the dynamics of external DM layers, as well. 

The most significant physical effect, however, is the outflow of baryon layers from the DM 
top-hat. The outflown baryon fraction increases with a. For a ~ 0.92, i.e. when DM and 
inner baryons have attained their virialization radius, the fraction of baryons which have 
leaked out from the fluctuation is so large as ~ 10 %, even for ft = 0.05; it reaches ~ 40 % 
for ft = 0.2 . These values increase by an additional 10 % in the RP case. 

Virialization 

In Figures El and ?? we assumed the present time to coincide with virialization for all 
DM inside the top-hat fluctuation and all baryons kept inside it, excluding therefore a large 
deal of baryons, either 10 or 40 %, initially inside the top-hat, but leaking out during the 
fluctuation growth. 

It must also be outlined that much care ought to be taken to define the virialization 
condition, for materials within any radius R, reading 



and taking into account that, during the whole fluctuation evolution, potential energies 
include three terms, due to self-interaction, mutual interaction, interaction with the DE 



2T{R) = RdU(R)/dR , 



(15) 



by summing up DM and baryon kinetic energies 




(16) 
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field. More in detail, for DM and baryons, we have 

U C (R) = U CC (R) + U cb (R) + U C>DE (R) = Art [* dr r 2 p c (r) [* c (r) + %(r) + * DE (r)\ , (17) 

Jo 

U b (R) = U bb (R) + U bc (R) + U b>DE {R) = 4vr F dr r 2 p b (r) [* 6 (r) + * c (r) + * DE (r)) , (18) 

*/ 

respectively. While 

Art Art 
Mr) = -^G Pb {r)r 2 , ^ DE (r) = -—Gp DE (r)r 2 , (19) 

in both expressions, a subtle difference exists for \l/ c . In U b we have simply 

Att 

*c(r) = - yG Pc (r)r 2 , (20) 

as for the other components, but this expression is different in U c , where it reads 

An 

Hr) = - y{7G[p c - p c {r)\ + Gp c }r 2 , (21) 

p c being the background DM density. The different dynamical effect of background and DM 
fluctuation arises from the different ways how its interaction with DE is treated. Energy 
exchanges between DM and DE, for the background, are accounted for by the r.h.s. terms 
in eqs. (jlj). In this case no newtonian approximation was possible and was made. For DM 
fluctuation, instead, the effects of DM-DE exchanges are described by a correction to the 
gravitational constant G, becoming G* = jG, which adds to the dependence of DM density 
on (p. This ought to be taken into account in the fluctuation evolution, as is done in eq. (|21j). 

The above expressions hold for any R. However, if R coincides with DM top-hat (or is 
smaller than it), the densities p b , c ,DE do not depend on r. Densities depend on r if one aims 
to take somehow into account the shells where baryons, initially belonging to the top-hat, 
have flown. In the former case, the integrals (|17|) . (|18j) and (|16j) can be easily performed 
and yield 

U C (R) = (3/5)G[M c (M b + M c ) + 1 M c ^M b }/R - (An /h)M c p DE R 2 (22) 

U b (R) = -(3/5)G[M b (M b + M c )}/R - (An/5)M b p DE R 2 (23) 
T C (R) = (3/10) M C R 2 (24) 
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where the relation r/r = R/R is used to calculate T C (R). Note that this relation is not 
valid for T b (R) because different baryon layers have different growth rates. Kinetic energy 
for baryons is then obtained by 

UR) = Y.n = Y.\K{K? (25) 

n n " 

Here the sum is extended on all R% < R. 

If integrals are extended above the DM top-hat, recourse to numerical computations is 
needed. 

MASS FUNCTIONS IN cDE THEORIES 

The expected physical behavior of a top-hat fluctuation, including non-linear features, 
can be compared with the evolution of that fluctuation if we assume that linear equation 
however hold, indipendently of the actual amplitude the fluctuation has reached. Actual 
gravitation prescribes that fluctuations approaching unity abandon the linear regime, slow 
down their expansion rate, reach maximum expansion, turn-around and recontract, finally 
recollapsing to nil. While this occur, we can formally assume that linear equation still hold 
and seek the value 5 rc that linear fluctuations would have at the time r rc when, according 
to actual gravitation, they have recollapsed. 

Let us also remind that the linear evolution does not affect fluctuation amplitude dis- 
tributions. Accordingly, if fluctuation amplitudes are distributed in a Gaussian way, this 
does not change in the linear regime. Therefore, at the time r rc , we can integrate on the 
distribution of fluctuation amplitudes for a given scale, taking those > 5 rc , so finding the 
probability that an object has formed and virialized over such scale. 

As already outlined in the previous section, full recollapse is not expected to occur. 
The usual assumption is however that the time running between the achievement of the 
virialization prescription and the time formally required for full recollapse is taken by the 
fluctuation to achieve a relaxed virial equilibrium configuration. This is the basic pattern 
of the PS-ST approach, that we aim now to apply to cDE cosmologies. In the presence 
of DM-DE coupling, however, a novel feature must be considered: starting from the initial 
amplitudes set by the linear theory, DM and baryon fluctuations require different times to 
reach full recollapse. Accordingly, if we require a baryon fluctuation to reach full recollapse 
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at the present time r (or at any time r), the initial linear 5b and 5 C must have been greater 
than those required to allow full recollapse at the present time r D (or at any time r) for a 
DM fluctuation. 

Let us also remind that the linear theory does not prescribe equal linear amplitudes 
for DM and baryons, but that their ratio is given by eq. ([lip. Therefore, it becomes also 
important to outline that we chose to refer to DM fluctuation amplitudes in the linear regime. 
Accordingly, let 5j$l c and 5^} c be both DM fluctuation amplitudes, defined so that, if the 
linear theory yields the former (latter) value at a time r, the corresponding DM (baryon) 
fluctuation has fully recollapsed at r. 
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FIG. 6: Ratio between DM and baryon mass in virialized halos, against the value used for S c ^ rc 
(DM masses are rescaled to the value they reach at r Q . The dashed line is the background Q, c /£lb 
ratio. 

In Figure |U we plot the (3 dependence of S^ c } c and 8fl c , at z = 0. They start from equal 
values for (5 = and gradually split as (3 increases. In Figure El we plot the z dependence of 
5(% and 6® c for (3 = 0.20 and 0.05. These values were computed by starting at z = 1000, 
taking into account also the radiative component, and using the full set of equations. 

Following the PS approach, the differential mass function then reads 

■m) = 2^ r ds M A- ( ' ««) 

M JSC,™ dM da M { v27rcr A / J 

d VM d ^ VM e-^. (26) 

V 71 M J5 c ,rc/CTM dM 

Here vu = 8m /&m an d 8 C ^ TC can be either 8^} c or 8f} c (or any intermediate value) according 
to Fig. 13] and 03 a choice shall be based on the observable to be fitted. 
The ST expression is obtainable from eq. ()26|) through the replacement 

v M exp(-i&/2) -> jV'v' m (1 + v'm /5 ) exp(-i/^/2) , (27) 
with JV' = 0.322 , v' M = 0.707 u 2 M , 
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meant to take into account the effects of non-sphericity in the halo growth. 

In the absence of coupling and baryon-DM segregation, the mass M in the PS and ST 
expressions (|26|) - (|27|) is the mass originally in the top-hat, which will then be comprised 
within a virial radius R v , such that M/(A7c/3)pRl = A v . 

In the presence of coupling and segregation the situation is more complex. However, in- 
dipendently of the value taken for <5 CjCr , the resulting virialized system will be baryon depleted. 
Different possible <5 CiCr 's will correspond to different depletions, but the final system shall 
however contain a smaller fraction of baryons, in respect the background VL b /VL c ratio. 

It is important to distinguish between two effects: (i) DM mass variation, (ii) The 
dynamics of gravitational growth. Let Mj and M vir be the masses, at the initial time and 
at virialization, rescaled to the values they will have at r G , so that the (i) effect is isolated. 
Then, while 

U m? • u\ (njnjMt + (jynjM, , (28) 

so that Mf/M b = fl c /flb , in the decomposition 

M mr = M c mr + M b vir (29) 

it will however be M c vir jM h mr > Q c /Q b . 

Let us now duly take into account also the (i) effect and consider the case when the 
mass function is set by 8^} c . Then, while M° ir /M? = exp[— C(4> V i r — 0j)], it will obviously 
be M b ir < M\: several baryon layers, initially belonging to the fluctuation, have not yet 
recollapsed or virialized. 

Let us then consider the mass function set by 5f} c . In this case it is M b vir = Mf, but it 
will be M^ ir /M![ > exp[— C(<p V i r — (pi)}. The extra DM mass is due to those layers, initially 
external to the fluctuation, first compressed and then conveyed inside the virialization radius, 
together with the baryons previously outflown from the DM bulk. 

For any <5 c rc in the S^ c } c -S^ b } c interval, some baryon layers will still be out and some extra 
DM will have been conveyed inside the virial radius by the fall out of outflown baryons. 
Hence, it will however be M£ ir exp[— C((f> — 4>vir] /M b ir > VL C /Vt b . 

In Figure El we plot this ratio, as a function of <5 Circ , in the S^ c -5^ b } c interval. The plot 
shows that, after a fast decrease, the ratio tends to a steady value, however exceeding the 
background ratio. The curves shown in this plot depend on the assumed (top-hat) shape 
for the primeval fluctuation, but similar curves would hold for any initial shape. 
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A prediction of cDE theories, therefore, is that Q c /Qb, measured in any virialized struc- 
ture, exceeds the background ratio. The excess is greater for structures where only DM has 
virialized. They might be characterized by an apparent disorder in the baryon component, 
still unsettled in virial equilibrium while, e.g., a lensing analysis would show that they are 
safely bound systems. 

We shall now plot mass functions obtained using either 5^} c or 5f} c . We expect ac- 
tual measures to yield a value comprised in this interval and, however, closer to the 6j?} c 
curve when baryon stripping is stronger. In Figure we plot the integral mass functions 
n(>M) = dM' tp(M') obtained through ST expressions (}27j l . Let us remind that 
large differences between models were never found in mass functions at z = 0, because of 
DE nature. The upper panel of each figure shows the mass function in the usual fashion, 
as often plotted to fit data or simulations. In the lower panel we plot the ratio between 
expected halo numbers for each model and ACDM. This confirms the small shifts between 
ACDM and dDE cosmologies, yielding just a slight excess, ~ 10 %, on the very large cluster 
scale, where observed clusters few units. 

Discrepancies can be more relevant between ACDM and cDE, whose effective mass func- 
tion should however lie inside the dashed areas, limited by the function obtained by inte- 
grating from 8^ c } c or 5fl c . The plots show a shortage of larger clusters. For (3 = 0.20, they 
are half of ACDM at ~ 3 ■ lO 14 ^" 1 !/© and really just a few above some lO 15 /i -1 M . Any 
realistic mass function, laying in the dashed area, can be falsified by samples just slightly 
richer than those now available. 

For (3 = 0.05 the shift is smaller, hardly reaching 20 %, still in the direction opposite to 
dDE. Here we meet what appears to be a widespread feature of cDE models: the discrepancy 
of dDE from ACDM is partially or totally erased even by a fairly small DM-DE coupling, 
and many cDE predictions lay on the opposite side of ACDM, in respect to dDE. If a ACDM 
model is then used to fit galaxy or cluster data arising in a cDE cosmology, we expect cluster 
data to yield a best-fit o~g value smaller than the one obtained by fitting galaxy data. 

The integral mass functions obtained through the original PS expression are shown in 
Fig. |HJ This figure shows just slight quantitative shifts. In the cases illustrated by the next 
figures, results from PS expression will be therefore omitted. 

Figure M compares RP results with SUGRA, at z = 0. Quantitative differences exist, 
when the self-interaction potential is changed, but most physical aspects are the same. 
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FIG. 7: Cluster per (Mpc /i -1 ) 3 , above mass M at z = 0, obtained from ST expression. Four 
models are considered, with equal f2's and h: ACDM, uncoupled SUGRA with A = 100 MeV, and 
two cDE with different /3's. The dashed areas are limited by the mass functions worked out for 
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FIG. 9: A comparison of SUGRA and RP results at z = using the ST expression. The lower 
panel of Fig. [7] is reproduced, neglecting dDE and = 0.05 cases and adding RP results (heavy 
dashed). 
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FIG. 10: Cluster per (Mpc /i -1 ) 3 , above mass M at different z's, in a fixed comoving volume, 
obtained from ST expression. Thick solid and dotted lines refer to cDE models. The dashed 
(thinner solid) line refers to dDE (ACDM). 
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FIG. 11: Evolution of the DM/baryon background ratio, due to the dynamics of the (ft field. 

In Figure El the redshift dependence of the expected cluster numbers in a comoving 
volume is plotted against the redshift z, for M = 10 14 and 4 • 1O 14 /i _1 M . As usual, the 
mass considered is the total cluster mass. As already widely outlined, the DM/baryon ratio 
in these masses however exceeds the background f2f,/f2 c ratio; the spread of the function 
corresponds to the spread of possible baryon/DM ratios. In top of that, however, one 
must also remind that the very background fi c /f^ ratio varies with redshift, because of the 
evolution of the <fi field. Hence, clusters observed at high z, in average, shall be however 
baryon poorer than present time clusters. The z dependence of the background fl c /^b ratio, 
for the model considered in this work, is plotted in Figure ITT1 

Plots similar to Figure El are often used to assert the possibility to discriminate between 
models and this plot is however significant to compare cDE with former results for other 
cosmologies. According to [17j, however, the discriminatory capacity of this observable can 
only be tested by plotting cluster numbers per solid angle and redshift interval, which also 
includes geometrical effects, often partially erasing dynamical effects. 

In the upper panels of Figures El and El cluster numbers per solid angle and redshift 
interval are plotted. In the lower panels we plot the ratios between each SUGRA model and 
ACDM mass functions. We consider again the mass scales 10 14 and 4 • lO 14 /i _1 M ; for the 
latter mass we provide a magnified box to follow the expected low-z behavior. Notice, in 
particular, that the high-;? behaviors for (3 = 0.05 or 0.2, for these mass scales, lay on the 
opposite sides of ACDM. The box in Fig. El show how we pass from numbers smaller than 
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FIG. 12: Number of clusters with M > 1O 14 /i _1 M in a fixed solid angle and redshift interval. 
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FIG. 13: Number of clusters above 4 • 10 14 /i 1 M & . Lines and comments as for Fig. 1121 
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FIG. 14: A comparison between SUGRA and RP results. Here the top panels of Figs. Upl and 1131 
are reproduced, omitting the j3 = 0.05 case and showing both SUGRA (dashed) and RP (heavy 
dashed) number intervals. The number distribution of ACDM is also given (solid line). 

ACDM to greater numbers, for (5 = 0.05, at a redshift z ~ 0.7. 

As is obvious, a better discrimination is attained for high masses or deep redshifts. How- 
ever, according to Fig. a sample including a few dozens clusters of mass > 4- lO 14 /i _1 M 
at z > 0.5 would already bear a significant discriminatory power. 

Altogether we see that, (i) when passing from ACDM to uncoupled SUGRA, the high-2: 
cluster number is expected to be greater, (ii) When coupling is added, the cluster number 
excess is reduced and the ACDM behavior is reapproached. (iii) A coupling (3 = 0.05 may 
still yield result on the upper side of ACDM, while (3 = 0.2 displaces the expected behavior 
well below ACDM. The ACDM behavior is approsimately met for (3 = 0.1. 

These behaviors arise, first of all, from the different evolution of pde at high z. In ACDM 
models, pbe simply keeps constant, while p m oc (1 + z) 3 , and therefore its relevance rapidly 
fades. In dDE models, instead, pde increases with z, although at a smaller rate than p m . 
To obtain the same amount of clusters at z = requires that they have existed since earlier. 

When coupling is added, however, gravitation is boosted by the $ field. In the newtonian 
language, this translates into a greater gravity constant (G*) and growing masses for DM 
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particles. This speeds up cluster formation and, when (3 increases, less clusters are needed 
at high z, to meet their present numbers. 

Finally, in Figure El we provide a comparison between the PS-ST mass function obtain- 
able for RP and SUGRA potentials. Only the case j3 = 0.2 is considered. A RP potential 
slightly strengthens the effects already seen for SUGRA; as a matter of fact, however, in 
most cases we find just small quantitative shifts. 

CONCLUSIONS 

In this work we aimed at predicting the cluster mass function in cDE cosmologies, by 
using the solution of the equations ruling the spherical growth of top-hat fluctuations in ST 
(or PS) expressions. The effectiveness of ST expressions has been widely verified for SCDM, 
ACDM and OCDM models. Also in simulations of models with dynamical (uncoupled) DE 
[loj l , ST expressions provide a fair fit of numerical outputs. 

The main finding of this study is the significant baryon-DM segregation, which has 
multiple effects. If cluster numbers are measured from their gravitational effects, e.g., by 
using lensing data, a PS-ST approach yields simple predictions. When cluster numbers are 
measured through other observables, we predict a number range, inside which observations 
should lie. The actual amount of objects, inside these ranges, is determined by a number of 
effects, that a PS-ST approximation cannot describe. 

An important example of such effects is the possible stripping of outer layers, in close 
encounters, which will mostly act on the baryon component. For a rather small coupling as 
(3 = 0.2, up to 40 % of the baryons belonging to the initial fluctuation could be stripped in 
this way and, even for the tiny coupling set by (3 = 0.05, 10% of baryons could be easily 
stripped. 

If cluster data are obtained from galaxy counts or hot gas features, they will exhibit the 
residual baryon amount. Indipendently of the baryon loss, which depends on the individual 
cluster history, cDE theories predict that the background fl c /^b ratio however increases with 
redshift. However, in top of that, the DM/baryon ratio measured in any virialized structure, 
exceeds the background ratio at the redshift where it is observed and is expected to exhibit 
significant variations in different systems, being smaller in larger systems, in average. 

It must be however outlined that the final baryon/DM ratio, in any galaxy cluster at 
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any redshift, even in the absence of any stripping effect, is expected to by smaller than the 
background fi&/f2 c ratio at that redshift. In this paper dedicated to a technical analysis 
of cDE mass functions we refrain from discussing this feature in further detail, although 
relating it with the apparent baryon shortage in clusters, (see, e.g., |18j) seems suggestive. 




Furthermore, when cluster data are obtained through the hot gas behavior, a complex 
interplay between baryons and potential well is expected. Once again, however, the model 
used for PS-ST estimates seems unsuitable to provide quantitative predictions, but anoma- 
lies in the temperature-luminosity relations are expected. 

A PS-ST analysis allows however to formulate further predictions, besides those con- 
cerning the flb/^c ratio in galaxy clusters and in galaxies. They concern the cluster mass 
function and its evolution. 

No large differences between models were ever found in the mass functions at z — 0, 
because of DE nature: just a slight excess, ~ 10 %, on the very large cluster scale, where 
observed clusters are a few units, was found in dDE, in respect to a ACDM cosmology. 

Discrepancies can be more relevant between ACDM and cDE, where a shortage of larger 
clusters is predicted. For (3 = 0.20, they are half of ACDM at ~ 3 • 10 u h~ 1 M Q and less than 
20 % above a few 1O 15 /i _1 M . Such strong shortage can be falsified by samples just slightly 
richer than those now available. For (3 = 0.05 the shift is smaller, hardly reaching 20 %, but 
still in the direction opposite to dDE. 

This is a widespread feature of cDE models: the discrepancy of dDE from ACDM is 
partially or totally erased even by a fairly small DM-DE coupling, and many cDE predictions 
lay on the opposite side of ACDM, in respect to dDE. Therefore, if a ACDM model (or any 
uncoupled model) is used to fit galaxy or cluster data arising in a cDE cosmology, we expect 
that cluster data may yield a smaller as, in comparison to the one worked out from other 
data sets. 

Turning to the evolutionary predictions, we expect an evolution faster than ACDM for 
any coupling (3 > 0.1, again the opposite of what we expect in uncoupled dDE. 

Most quantitative results given in this paper are worked out by assuming that the scalar 
field self-interaction potential is SUGRA, with A = 100 GeV. For the sake of comparison, 
in a few plots, results obtained for a RP potential are also shown. It should be reminded 
that, while the former potential predicts a linear behavior consistent with observations, the 
latter one can be fitted with linear observables only for quite low values of A, much below 
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the one considered here; we selected it just to provide a direct comparison tool, so allowing 
us to conclude that a different self-interacion potential may cause quantitative shifts up to 
some 10 %, but hardly affect our general conclusions. 

Altogether, there can be no doubt that cDE cosmologies open new prespectives for the 
solution of those problems where baryon-DM segregation due to hydrodynamics is appar- 
ently insufficient to explain observed features. These problems may range from the shortage 
of galaxy satellites in the local group, to galactic disk formation, up to the L vs T relation 
in large clusters. Stating how coupling can affect these and similar questions, by using a 
PS-ST approach, is hard. This calls for detailed n-body and hydrodynamical simulations 
of cDE cosmologies. 
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Appendix 1. The newtonian regime 



In the presence of inhomogeneities, the metric can read 

ds 2 = a 2 (r)[-(l + 2ip)dr 2 + (1 - 2ip)dxidx% (30) 

provided that no anisotropic stresses are considered, ip being the gravitational potential in 
the Newtonian gauge. Let us describe DE field fluctuations 5(p through 

<p = (4tt/3) 1/2 (5<j>/m p ) (31) 

and expand fluctuations in components of wavenumber k; let also be A = 71/ k. Let then be 

/ = 0-y3/167rG HV/V ), / 1 = ^| + /, / a = ^ + 2/ + / i; (32) 

V Q being a reference value of the potential. It is also useful to define Y 2 = 8nGV((f)) a 2 /3H 2 . 
The equations ruling the evolution of the ip field and gravity, keeping just the lowest order 
terms in A, as is needed to obtain their Newtonian limit, then read 

ifj = --\ 2 ({l b 5 b + tt c 5 c + 6X<p + 2Xtp' -2Y 2 f\<p) , ip' = 3x<p - ip , (33) 

h' 

y" + (2 + —)<f? + \- 2 <p - 12X^ + 4i(jX + 2Y 2 (f 2 ^ - M) = (3Q C (5 C + 2iP) ; (34) 
rt 

let us remind that X is defined in eq. (jSJ) and notice that, if DE kinetic (and/or potential) 
energy substantially contributes to the expansion source, X (and/or Y) is 0(1). 

In the Newtonian limit, <p derivatives shall be neglected, the oscillations of ip and the 
potential term f 2 Y 2 (p should be averaged out, by requiring that A << (f 2 Y)~ 1 , and, in 
eq. (JBH), the metric potential ip (ex A 2 ) can also be neglected. Then, eqs. (jSSj) ~ (iHl|) become 

ip = ~\\n h 5 b + n c s c ), , rV«. ■ (35) 

Baryon and DM density fluctuations are then ruled by the eqs. derived from the 

stress-energy pseudo-conservation T^.^ = 0. Thereinside, ' yields differentiation with respect 
to In a. [Let us then notice that, taking << fi c , putting 5 C oc e-> > x ^> da and 5b = b5 c with 
6 = cost, eqs. (jBJ), allow us to obtain the bias factor fTTf ]. 

The acceleration of a DM or baryon particle of mass m c ^ can then be derived from 
eqs. ©-(0). Let it be in the void, at a distance r from the origin, where a DM (or baryon) 
particle of mass M c (or M b ) is set, and let us remind that, while p b oc a -3 , it is 

Pc = p oc a- 3 e" c ^-^), p Mc = M oc a~ 3 e- C ^-^5(0), (36) 

26 



because of the DE-DM coupling (the subscript o indicates values at the present time r G ; 
let also be a Q = 1). We can then assign to each DM particle a varying mass M c (0) = 
M oc e~ c ^-M 

Then, owing to eq. (|36j) . and assuming that the density of the particle widely exceeds the 
background density, it is 

V C 8 C = = |^M C (0)5(O) , n b 5 b = = |^M^(0), (37) 

per 6ri A a p cr 6n z a 

(p cr is the critical density and S is the Dirac distribution). Reminding that V • v Cj & = 9 cb H 
and using the ordinary (not conformal) time, eq. © yields 

V ■ v c = -H(l - 2f3X) V ■ v c - ArrGa- 2 {^M c {<j>) + M b ) 5(0) (38) 

(dots yield differentiation in respect to ordinary time and H = a /a). Taking into account 
that the acceleration is radial, as the attracting particles are in the origin, it will be 



J d 3 r V ■ v = 47r J dr d(r 2 v)/dr = A-Kr 2 v. 



Accordingly, the radial acceleration of a DM particle is 



n ory\u G*M c (<j)) GM b 
v c = -(1 - 2f3X)Hv c ■ n — — , (39) 

(n is a unit vector in the radial direction; R = ar). 

By following a similar procedure for baryons, we obtain: 

GMM GM b 
v b = -HY b ■ n — — (40) 

In the presence of full spherical symmetry, Vb, c • n = v b>c \ being then b = r b and c = r c , 
eqs. (|14|) immediately follow. 

For the sake of completeness, here below, we write also the equations for the physical 
DM and baryons radii which easily follow from eqs. ()39j) and ()40j) . While for baryons we 
have the usual Friedman-like equations 

Air 

Rb = ~y g + P b + P DE ^ + 3w )] R b ' ( 41 ) 

for DM we have 

R n c = C4>R n c - C<j)HR n c - y G [p c + 7 (p c - Pc) +Pb + Pde(1 + 3u;)] R n c . (42) 
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